NSWC/WOL/TR  75100 


O 

<0  * 
(A 


\JFCHMICAL 

r^REPOR^T1 

uu 


WHITE  OAK  LABORATORY 


CQ  PREDICTION  OF  VERTICAL  WATER-ENTRY  FORCES  ON  OGIVES  FROM  CONE  DATA 


BY 

^ J.L.  Baldwin 


NAVAL  SURFACE  WEAPONS  CENTER 
WHITE  OAK  LABORATORY 
SILVER  SPRING,  MARYLAND  20910 


• Approved  for  public  release;  distribution  unlimited. 


15  AUGUST  1975 


DOC  . 

\\~  1Q7«  1 ! i 

Jj  APR  *3  ,3?0  ' ; 0 

...  • ill; 

1 UlafcJfc-'-  u 


NAVAL  SURFACE  WEAPONS  CENTER 

WHITE  OAK,  SILVER  SPRING,  MARYLAND  20910 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  of  THIS  PAGE  (When  Date  Entered) 


n 

n 


KSi 

(i 


p 


arr 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


■ — 1 REPORT  |NUMBUER^*WOTfwwu^i^‘*'f',y'v\ 

t/.  ) NSWC/w6l/TR«»7  5-1^0/ 


2.  GOVT  ACCESSION  NO.I  3,  RECIPIENT’S  CATALOG  NUMBER 


Prediction  of  Vertical  Water-Entry 
Forces  on  Ogives  from  Cone  Data# 


&I..tuSEM«e*€OVER,ED 


^nterim 


«|.^tRrORMfnff'6T?G:-rtfP6RTNU 


B.  CONTRACT  OR  GRANT  NUMBERf*; 


J.  L. /Baldwin 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 


ECT,  TASK 


— < — AR*A-«r-woRtrtmfT*NimtjefA 

Naval  Surface  Weapons  Center  / 

White  Oak  Laboratory  ^ l*y|  _ — ~4|- 

White  Oak.  Silver  Spring.  Maryland  20910 


II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 


pnaT-nATC. 


15  Aug 


T«  MONITORING  AGENCY  NAME  A ADDRESS  (It  dltlerent  from  Controlling  Ottlce) 


MBtrrarvkfSis" 

48 


ILL 


15.  SECURITY  CLASS,  (ot  thla  report) 


-line  lass  ifiad 

. DECL  ASSIFrCATlON/TlC 


t5*.  DECLASSlFrCATlONADOWNGRADING 
SCHEDULE 


16.  DISTRIBUTION  STATEMENT  (ol  Me  Report) 

Approved  for  public  release;  distribution  unlimited 


<7.  DISTRIBUTION  STATEMENT  (ol  the  abatrecl  entered  In  Block  20,  II  dlllarent  from  Report) 


IB.  SUPPLEMENTARY  NOTES 


19.  KEY  WOROS  (Continue  on  ravarae  aide  II  necaaaery  end  Identity  by  block  number) 

Ogives,  Cusps,  Cones,  Water  entry,  Vertical  drag 


w 


ABSTRACT  (Continue  on  ravarae  aide  It  necaaaery  and  Identity  by  block  number) 

An  empirical  method  of  predicting  the  drag  coefficient  for 
vertical  water  entry  of  pointed  shapes  is  presented.  The 
drag  coefficient  versus  depth  function  is  generated  from 
experimental  results  of  the  vertical  entry  of  cones  of  various 
angles.  Comparison  with  experimental  results  shows  that  the 
predicted  values  are  accurate  to  within (±2 5 percent.  A computer 
program  for  predicting  the  water-entry  djqag_x:oeff icient  of  many 
shapes  is  included,  a 


DD  i jan^73  1473 


EDITION  OF  I NOV  SB  IS  OBSOLETE 
S/N  0 102*014*  66o\  | 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  Dele  tntered) 


S 


CONTENTS 


Page 

INTRODUCTION 5 

BACKGROUND 7 

PREDICTION  METHOD 7 

COMPARISON  WITH  EXPERIMENTS 9 

SUMMARY  AND  CONCLUSIONS 15 

REFERENCES 16 

APPENDIX  A A-l 

APPENDIX  B B-I 

APPENDIX  C C-.L 

ILLUSTRATIONS 

Figure  Title  Page 

1 Expanding  Disk  Predictions  for  90  and  20- 
Degree  Cones  Compared  to  Experimental 

Results 18 

2 Water-Entry  Nomenclature 19 

3 Equivalent  Cones 70 

4 Total  Drag  Coefficient  Vs  Depth  for  60/43 

Ogives 21 

5 Total  Drag  Coefficient  Vs  Depth  for  88.4/ 

32  Ogives 22 

6 Total  Drag  Coefficient  Vs  Depth  for  59.94/ 

-43.4  Cusps 23 

7 Total  Drag  Coefficient  Vs  Depth  for  88.94/ 

-31 24 

8 Total  Drag  Coefficient  Vs  Cone  Angle  for 

Tangent  Ogives 25 

9 Total  Drag  Coefficient  Vs  Swept  Angle  for 

6 0/B  Ogives 26 

A-l  Ogive  Geometry A-4 

A-2  Base  Area  Comparison  of  Secant  and  Tangent 

Ogives A-5 

B-l  Listing  and  Sample  Run  of  Computer  Program 

Written  in  BASIC B-4 

C-l  Total  Drag  Coefficient  Vs  Depth  for  90/90 

Ogives  (Hemispheres) C-2 

C-2  Total  Drag  Coefficient  Vs  Depth  for  60/60 

Ogives C-3 

C-3  Total  Drag  Coefficient  Vs  Depth  for  48.18/ 

48.18  Ogives C-4 

C-4  Total  Drag  Coefficient  Vs  Depth  for  41.44/ 

41.44  Ogives C-5 


2 


NSWC/WOL/TR  75-100 


LIST  OF  SYMBOLS 


Base  area  (ttR2) 
Buoyancy  force 


A constant 


Total  drag  coefficient 
Steady-state  pressure  drag  coefficient 
Transient  drag  coefficient 
Friction  drag  coefficient 


Missile  diameter 


Gravitational  constant 


Penetration  ratio 


rS  at  maximum 
i cone  length. 


Total  mass  constant 


Added  mass  constant 
Length 

Missile  mass 


Added  mass 


Base  radius 


Local  radius 


Distance  traveled  after  water  contact 


Instantaneous  speed 
Water  contact  speed 
Total  cone  angle 
Swept  angle 
Fluid  mass  density 


NSWC/WOL/TR  75-100 
INTRODUCTION 


The  two-dimensional  case  of  water  impact  originally 
received  the  majority  of  attention  due  to  its  application 
to  the  landing  of  flying  boats  and  seaplanes.  Analytical 
prediction  of  the  forces  resulting  from  water  entry  began 
with  a study  of  the  landing  impact  of  seaplane  floats  by 
von  Karrnan  in  1929  (Ref.  (1)).  In  his  paper,  von  Karman 
assumed  that  all  the  impact  force  was  due  to  an  increase 
in  added  mass,  and  that  the  amount  of  added  mass  was 
identical  to  one  half  of  the  added  mass  of  a flat  strip 
with  width  equal  to  the  width  of  the  float  where  inter- 
sected by  the  original  water  surface.  The  factor  of  one 
half  was  used  because  half  of  the  assumed  flow  field  was 
in  air  which  contributed  a negligible  amount  of  mass  to 
the  total.  In  1931,  Wagner  considered  the  effects  of 
water-surface  displacement  during  entry  (Ref.  (2)).  Many 
subsequent  investigators  have  also  considered  the  problem 
of  surface  motion.  Mayo  in  1945  (Ref.  (3))  and  Pierson 
in  1951  (Ref.  (4))  were  still  refining  the  prediction  of 
seaplane  impact.  In  addition  to  their  specific  contri- 
butions, they  also  include  summaries  of  much  of  the  previous 
work  on  the  seaplane  problem. 


The  prediction  of  the  water-impact  force  experienced  by 
axisymmetric  bodies  began  as  early  as  1942  when  Plesset 
applied  von  Karman' s approach  to  the  entry  of  torpedoes. 

In  this  method,  the  assumption  was  made  that  the  added 
mass,  "m"  at  any  point  during  entry  was  the  same  as  half  that 
of  a thin  flat  circular  plate  defined  by  the  intersection 
of  the  nose  and  the  original  water  surface.  Half  of  the 
added  mass  for  a disk  is 

m = 1.28  p r3  (1) 


where  r is  the  radius  and  p is  the  water  density.  Differentiation 
with  respect  to  the  distance  s yields: 

= (1.28)  p irr2  % (2) 

ds  as 

Also,  the  transient  drag  coefficient  "C^t"  had  been  shown  to  be 
given  by  equation  (3)  (See  Appendix  A). 


r = _i  dm 
pA  ds 


(3) 


Substituting  equation  (2)  into  equation  (3)  yields: 
cdt 

where  R is  the  base  radius  and  A is  the  base  area. 


2J>£  vr2  dr  = 256 
A ds 


(r)2  dr 
\R/  ds 


(4) 
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For  example;  for  a cone  of  total  angle  "a": 

r = S tan  (j'j  ; Cdt  =2.56  ^tan  ^ for  0 < S < R tan  f 

Cj.  - 0 for  S > tan  R — 
at  2 

For  a nose  shape  following  a cubic  relationship 

r = CS1//3;  Cdt  = — 8-5-C-3  for  0 < S < (R/C)3 

r2 

Cdt  = C for  S > (R/C)3 

Then  Cdt(S)  was  easily  calculated  for  any  desired  ^hape;  however, 
as  the  theory  was  developed  and  experimental  results  accumulated, 
several  substantial  deficiencies  became  apparent.  The  predicted 
drag  coefficients  for  90  and  20  degree  total  angle  cones  are  com- 
pared to  experimental  values  (Ref.  (5))  in  Figure  1.  Experi- 
mentally derived  water-entry  drag  coefficients  for  cones  are  also 
given  in  Reference  (6)). 

The  drag  force  predicted  by  this  method  becomes  equal  to 
zero  when  the  base  of  the  nose  passes  the  undisturbed  water 
surface  and  remains  at  zero  for  deeper  penetration.  In  reality, 
the  transient  axial  force  continues  below  the  water  surface  for 
about  a nose  length  before  the  steady-state  drag  coefficient  is 
obtained. 

In  most  cases,  the  predicted  drag  coefficient  is  much 
larger  than  that  obtained  experimentally.  This  over-prediction 
seems  to  be  caused  by  the  assumption  that  the  nose  shape  does 
not  affect  the  added  mass  constant.  Note  that  this  method 
predicts  that  the  final  added  mass  and  hence,  the  change  in 
velocity  during  water  entry  is  a function  of  only  missile 
diameter. 

Added  mass  constants  for  the  various  cones  were  derived 
theoretically  by  Shiftman  and  Spencer  (Ref.  (7))  and  experi- 
mentally (Ref.  (5)).  In  both  cases,  the  mass  constant  was 
equal  to  zero  at  zero  cone  angle  and  increased  monotonically 
with  increasing  cone  angle. 

The  predicted  value  of  rise  time  is  normally  longer 
than  experimental  value.  This  error  results  from  the  fact 
that  the  water  surface  usually  does  not  remain  flat  during 
entry  as  assumed. 
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It  is  the  purpose  of  this  report  to  present  a simple 
systematic  method  of  applying  the  experimental  results 
for  cones  given  in  reference  (5)  to  the  prediction  of  the 
total  drag  coefficient  experienced  during  vertical  water 
entry  by  ogives  and  other  non-flat  shapes. 

BACKGROUND 

The  basic  equations  for  axial  forces  encountered  during 
water  entry  are  included  in  Appendix  A,  which  was  copied  from 
reference  (5).  The  identification  of  ogives  in  this  report 
follows  the  method  used  in  reference  ( 8 ) , each  ogive  shape 
being  uniquely  defined  by  a cone  angle  and  a swept  angle. 

The  ogive  geometry  section  of  reference  ( 8 ) is  also  included 
in  Appendix  A. 


PREDICTION  METHOD 

The  proposed  technique  for  estimating  the  total  drag 
coefficient  of  axisymmetric  bodies  during  vertical  water 
entry  consists  of  predicting  a friction  drag  coefficient, 
a steady-state  pressure  drag  coefficient  and  a transient 
drag  coefficient,  and  then  summing  the  three  functions.  Or 
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Cd(S)  = Cf  (S)  + Cds(S)  + Cdt(S)  (5) 


where  each  coefficient  is  for  constant  velocity  entry  and 
constant  velocity  cone  entry  data  is  used  to  calculate  Cds 
and  Cdt«  A computer  program  written  in  BASIC  that  performs 
these  calculations  is  presented  in  Appendix  B. 

During  the  wetting  process  the  water-entry  force  acts 
upon  an  area  that  extends  above  the  original  water  surface 
as  shown  in  Figure  2.  The  upper  edge  of  this  area  is  defined 
as  the  effective  water  surface  and  its  distance  from  the  tip 
of  the  missile  is  defined  as  the  effective  distance  (Seff) . 
When  the  shape  is  a cone,  self  similar  flow  exists  which 
implies  that  the  distance  from  the  tip  to  the  original  water 
surface  divided  by  Seff  produces  a constant.  This  ratio  is 
approximated  by  the  penetration  ratio  (h) . 


S at  maximum  cd  _ , ^ S (6) 

Cone  Length  ” ~ Sef£ 

The  values  of  penetration  ratio  for  total  cone  angles  between 
ten  degrees  and  140  degrees  was  given  in  reference  (5)* 
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For  more  general  shapes  the  value  of  Seff  is  estimated  by 
adding  the  incremented  effect  of  each  body  angle  as  it  passes 
the  original  water  surface  as  shown  by 

seff  = f ^ (7) 

l ” 

where  h = f(a)  and  a = g(S).  Wetting  is  assumed  completed  when 
the  effective  distance  equals  the  length  of  the  nose. 

The  friction  drag  coefficient  is  estimated  by:  integrating 

a constant  friction  coefficient  of  .003  over  the  body  to  the 
effective  water  surface;  taking  the  axial  component  and  correcting 
the  reference  area  to  the  base  area  of  the  body.  This  simple 
form  was  chosen  so  that  the  total  drag  coefficient  would  not 
be  a function  of  entry  speed  and  size.  The  addition  of  a varying 
friction  coefficient  could  be  made  to  the  computer  program  sub- 
routine beginning  at  line  1200  for  increased  accuracy  for  those 
cases  that  have  significant  friction  drag. 

The  steady-state  pressure  drag  coefficient  is  estimated 
to  be  the  same  as  a cone  with  the  same  base  diameter  as  the 
missile,  undisturbed  water  surface  diameter,  and  the  same 
effective  water  surface  body  angle  as  the  missile  as  illus- 
trated in  Figure  3.  Addition  of  the  friction  drag  coefficient 
to  the  steady-state  pressure  drag  coefficient  produces  the 
steady- state  drag  coefficient  which  is  computed  up  to  the 
first  maximum  larger  than  .02.  At  all  deeper  depths  the 
maximum  value  is  used.  When  no  maximum  occurs  on  the  body, 
the  last  value  computed  during  wetting  is  defined  as  the 
maximum,  and  applied  as  above.  In  effect  cavitation  is  pre- 
dicted to  begin  at  the  depth  which  maximizes  the  steady-state 
drag  coefficient. 

The  transient  drag  coefficient  is  also  estimated  from  an 
equivalent  cone.  This  cone  has  the  same  base  diameter  as  the 
missile.  The  equivalent  cone  angle  is  the  larger  of  the  two 
body  angles  taken  at  the  undisturbed  water  surface  and  at  the 
effective  water  surface  (see  Figure  3) . When  wetting  is  com- 
pleted the  last  value  of  transient  drag  coefficient  and 
equivalent  cone  angle  are  retained  in  the  program.  At  deeper 
depths  the  transient  drag  coefficient  decays  at  the  rate  given 
in  ref erence  (5) for  the  equivalent  cone  angle.  When  the  value 
of  becomes  smaller  than  two  percent  of  C^g,  all  calculations 
are  terminated. 
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The  calculating  procedure  may  seem  to  be  arbitrary  and 
to  violate  accepted  theory;  however,  it  is  presented  as  an 
engineering  tool  and  during  its  development  this  method 
resulted  in  the  best  overall  fit  to  available  experimental 
data. 


COMPARISON  WITH  EXPERIMENTS 

The  predicted  values  of  the  drag  coefficient  vs.  depth 
function  were  computed  for  all  axisymmetric  shapes  for  which 
water-entry  data  were  found.  In  addition,  the  steady-state 
drag  coefficient  was  computed  for  all  shapes  for  which  the 
experimental  value  was  found.  Comparison  of  the  calculated 
values  with  the  experimental  data  showed  that  the  predicted 
values  would  satisfy  the  accuracy  requirements  of  most 
engineering  applications. 

The  results  of  reference  (5)  were  used  to  develop  the 
numerical  curve  fits  that  form  the  core  of  the  computing 
program.  In  order  to  determine  the  accuracy  of  the  fitting 
operation,  the  maximum  total  drag  coefficient,  the  penetration 
ratio  and  the  final  steady-state  drag  coefficient  were  com- 
puted for  each  cone  angle  reported  in  reference  (5)  and  also 
for  several  blunter  cones . The  results  of  these  calculations 
are  compared  to  the  experimental  results  in  Table  1. 

In  the  actual  calculation  procedure,  the  instantaneous 
drag  coefficients  are  used;  hence,  it  is  of  the  utmost 
importance  to  have  an  accurate  fit  with  regard  to  penetration 
ratio.  An  inspection  of  the  results  shows  that  the  penetration 
ratio  is  calculated  to  within  one  percent  of  the  experimental 
data  for  total  cone  angles  between  20  and  140  degrees.  Pene- 
tration ratios  larger  than  unity  were  not  allowed  in  the  fitting 
process,  which  accounts  for  the  large  errors  for  the  10“ and 
15-degree  cones.  These  large  errors  in  penetration  ratio  for 
the  ten  and  15-degree  cones  caused  the  total  maximum  drag 
coefficient  and  the  steady-state  drag  coefficient  to  also  be 
lower  than  expected.  In  general,  the  values  of  calculated 
maximum  total  drag  coefficients  were  larger  than  the  experi- 
mental values.  As  indicated  in  reference  (5) , this  is  to  be 
expected  due  to  "roll  off"  that  occurs  at  the  peaks. 

Upon  discarding  the  values  for  maximum  total  drag  coef- 
ficient for  30-degree  cones,  the  remaining  calculations  and 
experimental  values  agree  to  within  + ten  percent. 

Steady-state  drag  coefficient  values  were  found  for  ogives 
and  paraboloids  of  revolution.  The  comparison  of  the  experi- 
mental values  at  zero  cavitation  number  and  the  calculated 
values  are  given  in  Table  2 for  ogives  and  Table  3 for  paraboloids. 
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The  calculated  values  of  steady-state  drag  coefficients  for 
ogives  are  about  12  percent  higher  than  the  average  experi- 
mental values.  It  is  interesting  to  note  that  a cavitation 
number  change  of  o = .015  would  reduce  this  difference  by 
approximately  half.  Any  small  error  in  cavitation  number 
would  not  significantly  reduce  the  differences  for  the 
blunter  shapes. 

The  scant  experimental  data  available  pertaining  to  the 
steady-state  drag  coefficients  of  paraboloids  of  revolution 
are  compared  to  calculated  values  in  Table  3 . 

TABLE  3 

Comparison  of  Steady-State  Drag  Coefficient 
for  Paraboloids  of  Revolution 
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Experimental 

Calculated 

Cal./ 

L/D 

cds  (max) 
(Form) 

Ref. 

^ds (max) 
(Form) 

Exp. 

Form 

1 

.125 

15 

.138 

.130 

1.040 

3.33 

.025 

15 

.0415 

.0150 

.600 

Agreement  is  much  better  for  the  blunt  paraboloid  than  for 
any  of  the  blunt  ogives.  The  discrepancy  for  the  fine  para- 
boloid may  be  due  to  the  method  used  in  reference  (15)  to 
adjust  the  experimental  data  to  zero  cavitation  number.  How- 
ever, when  the  friction  drag  coefficient  is  added  to  the  form 
drag  coefficient  to  obtain  the  steady-state  drag  coefficient, 
the  difference  is  less  than  20  percent. 

The  principal  purpose  of  this  report  is  to  propose  a 
method  for  computing  the  total  drag  coefficient  for  pointed 
axisymmetric  shapes  during  vertical  water  entry.  In  order  to 
evaluate  the  accuracy  of  this  technique,  the  total  drag 
coefficient  versus  depth  of  penetration  function  was  com- 
puted for  all  relevent  shapes  for  which  water-entry  data  was 
available.  The  results  for  cones  have  been  previously  dis- 
cussed. The  calculated  values  for  ogives  and  cusps  are 
graphically  compared  with  the  experimental  data  in  Figures 
4-7  and  in  Appendix  C.  The  impact  data  for  other  shapes  such 
as  paraboloids  were  not  found.  However,  current ' studies  on 
such  shapes  should  yield  such  data  in  the  near  future. 
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Numerical  comparisons  were  easily  possible  with  the 
four  shapes  reported  in  reference  ( 8 ) . The  same  instru- 
mentation and  type  of  models  were  used  in  this  study  as 
were  used  in  the  study  of  cones,  which  gave  the  basic 
emperical  data  for  this  method.  Hence,  if  the  method  has 
any  usefulness,  one  would  expect  good  agreement  between  the 
calculated  functions  and  the  experimented  results.  The 
comparison  of  the  values  of  maximum  total  drag  coefficients 
and  depth  at  which  the  drag  coefficients  become  maximum  are 
shown  in  Table  4. 

It  is  seen  that  for  cusps  the  calculated  value  for  the 
maximum  drag  coefficient  is  about  12  percent  higher  than 
in  the  experimental  results . This  is  somewhat  greater  than 
the  error  in  reproducing  the  cone  data.  Perhaps  the  cusps 
have  a larger  percent  roll-off  than  cones.  The  ogives  are 
both  computed  about  four  percent  low,  which  is  interesting 
because  the  steady-state  drag  coefficients  for  these  shapes 
were  computed  higher  than  the  experimental  values.  In  each 
; of  tiie  four  cases,  the  computed  depth  at  which  the  maximum 

r drag  coefficient  occurs  agreed  with  the  experimental  data 

: to  within  four  percent.  Figures  4,  5,  6 and  7 show  the 

graphical  comparison  between  the  computed  drag  coefficients 
and  the  experimental  results  for  these  four  shapes. 

f 

; A comparison  of  computed  maximum  drag  coefficients  with 

| experimental  results  for  various  tangent  ogives  is  shown  on 

f Figure  8.  The  ogives  reported  in  reference  ( 8 ) were  all 

| secant  ogives  and  were  corrected  to  tangent  ogives  by  the 

| ratio  of  reference  areas.  As  with  the  numerical  comparison 

j shown  in  Table  4 , the  computed  values  of  the  maximum  total 

j drag  coefficients  seem  to  be  slightly  lower  than  most  of  the 

experimental  results.  The  sudden  increase  in  computed  drag 
coefficient  values  is  caused  by  friction.  This  effect  is 
accentuated  because  at  a total  cone  angle  of  about  12  degrees 
t and  finer,  the  steady-state  drag  coefficient  does  not  have 

1 a maximum  value  on  the  body.  Hence,  the  friction  drag  is 

applied  to  the  entire  surface  of  the  forebody.  This  result 
implies  that  the  forebody  becomes  fully  wetted  for  cone 
angles  of  12  degrees  and  less.  There  has  been  no  experi- 
mental verification  of  this  prediction  of  wetting;  therefore, 
the  predicted  drag  values  for  very  fine  ogives  should  be 
used  with  care. 

The  graph  of  calculated  maximum  total  drag  coef f,icients 
for  ogives  and  cusps  with  a 60-degree  cone  angle  is  shown  on 
, Figure  9.  Also  shown  are  the  experimental  values  for  four 

such  shapes.  As  before,  the  calculated  values  for  the  cusp 
and  cone  are  too  large,  while  the  calculated  values  for  the 
ogives  are  probably  too  small. 
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SUMMARY  AND  CONCLUSIONS 

An  empirical  method  of  predicting  the  drag  coefficient  for 
vertical  water  entry  of  pointed  shapes  was  developed.  This  method 
applied  the  experimental  results  previously  obtained  for  the 
vertical  water  entry  of  cones  of  various  angles  to  the  geometric 
shape  of  the  forebody  of  interest.  Predicted  drag  coefficient 
functions  agreed  to  within  ±25  percent  of  the  available  experimental 
results.  The  computer  program,  written  in  BASIC,  was  used  to 
calculate  the  predicted  drag  coefficient  function;  it  is  given 
in  Appendix  B.  For  those  who  wish  to  write  their  own  program, 
the  method  of  prediction  consists  of  the  following  steps: 

1.  Determine  an  expression  for  the  instantaneous  transitory 
drag  coefficient  for  all  cone  angles. 

2.  Determine  an  expression  for  the  instantaneous  forebody 
pressure  drag  coefficient  for  all  cone  angles. 

3.  Determine  an  expression  for  penetration  ratio  for  all 
cone  angles. 

4.  Define  the  zone  of  action  to  lie  between  the  original 
water  surface  and  the  effective  water  surface. 

5.  Divide  the  forebody  into  computing  stations  beginning 
at  the  tip. 

6.  For  each  station,  compute  the  effective  surface  position 
by  adding  the  effect  of  penetration  ratio  for  all  stations  to  the 
tip. 


7.  Calculate  the  friction  drag  from  the  tip  to  the  effective 
surface  associated  with  each  station,  and  convert  to  coefficient  form. 

8.  Calculate  the  pressure  drag  coefficient  based  upon  the 
angle  of  the  forebody  at  each  effective  surface  position. 

9.  Add  the  friction  drag  coefficient  and  the  forebody  pressure 
drag  coefficient  to  obtain  the  steady-state  drag  coefficient. 

10.  Use  the  value  of  the  first  maximum  of  the  steady-state  drag 
coefficient  for  all  stations  at  greater  depths. 

11.  At  each  station,  calculate  the  transitory  drag  coefficient, 
using  the  maximum  body  angle  in  the  zone  of  action. 

12.  Add  the  steady-state  drag  to  the  transitory  drag  coefficient 
to  obtain  the  total  drag  coefficient  at  each  station. 

13.  At  stations  for  which  the  forebody  base  has  been  wetted, 
reduce  the  transitory  drag  coefficient  proportionally  to  the  cone 
of  the  last  calculation  angle. 
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DEPTH/CONE  LENGTH 


FIG.  1 EXPANDING  DISK  PREDICTIONS  FOR  90  AND  20  DEGREE  CONES 
COMPARED  TO  EXPERIMENTAL  RESULTS 
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FIG.  4 TOTAL  DRAG  COEFFICIENT  VS  DEPTH  FOR  60/43  OGIVES 


EXPERIMENTAL  REF  (8  ) 
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FIG.  5 TOTAL  DRAG  COEFFICIENT  VS  DEPTH  FOR  88.4/32  OGIVES 


PREDICTED 


FIG.  6 TOTAL  DRAG  COEFFICIENT  VS  DEPTH  FOR  59.94/-43.4  CUSPS 


EXPERIMENTAL  REF  (8  ) 
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FIG.  8 TOTAL  DRAG  COEFFICIENT  VS  CONE  ANGLE  FOR 
TANGENT  OGIVES 
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FIG.  9 TOTAL  DRAG  COEFFICIENT  VS  SWEPT  ANGLE  FOR  60/e  OGIVES 


APPENDIX  A 


BASIC  EQUATION 


The  equation  relating  the  momentum  just  before  first 
contact  with  the  momentum  at  some  later  time  is  given  as 
Equation  (1) , and  the  force  equation  is  given  as  Equation  (2) 


U0M-U (M+m)  = C B dt-Mgt+P  C CdsAU2dt 


- ||  (M+m)  - U |!  = B-Mg  + t Cds  AU: 


um 

Substituting  U2  ds  for  the  second  term  of  Equation  (2)  and 
rearranging  terms  gives: 


||  (M+m)  - B + Mg  = ||  U2  + f Cds  AU2 


If  the  total  drag  coefficient  is  defined  as 


_ 2 dm 

C = + c 

d pA  dS  ds 


then  Equation  (3)  may  be  rewritten  as 


dU  1 2 

dt  (M+m)  - B + Mg  = - p Cd  AUZ 


The  value  of  added  mass  (m)  and  dS  can  be  calculated  using 
Equation  (3)  and  the  total  drag  coefficient  Cd  was  computed 
from  Equation  (4) . However,  a simpler  procedure  is  to 
determine  m from  (1)  and  from  (5) . In  order  to  ensure 
an  accurate  evaluation  of  the  integrals  in  Equation  (1) , 
the  added  mass  (m)  was  also  computed  by  integrating  dm  as 
given  in  Equation  (4)  using  the  computed  value  of  Cd9® 

For  the  vertical  entry  of  cones,  the  added  mass  should 
be  given  by 


jn  = pkr' 


pk (tan  3 3s3 
2 


Then 


dm  ■>  , ,*.  mx  3 _2 

g^r  = 3pk  (tan  J?)  S* 

2 


Substituting  (7)  into  (3)  yields 


UU  i y r\  * 

- j£  (M+m)  - B + Mg  * (3pk(tan  &)*S  + P cdsA)U‘ 

2 2 


If  the  total  mass  constant  (K)  is  defined  as: 

cds 

K = k + “g~*{tan  Q)  7r 

2 

then  (8)  can  be  rewritten  as 

dU  m 3 o 2 

- (M+m)  ~ B + Mg  = 3pK (tan  S)  s U 


This  form  of  the  equation  was  chosen  to  compare  easily 
the  experimental  results  with  the  theoretical  results 
obtained  in  reference  3. 

It  should  be  noted  that  the  total  mass  constant  (K) 
is  also  a drag  coefficient  related  to  the  area  intersected 
by  the  original  water  surface.  The  relationship  between 
the  mass  constant  (K)  and  the  total  drag  coefficient 
is  obtained  by  subtracting  Equation  (5)  from  Equation  (10) 
and  dividing  out  p and  U2: 

Cg  = 6K(tan  ^ (lj 

2 


If  the  area  used  to  relate  Cd  to  drag  force  is  the  area 
intersected  by  the  heaved  water  surface,  then  the  instan- 
taneous drag  coefficient  is  given  by 


'd^  = (6K  tan  ^-)  (ir) 


(12) 


OGIVE  GEOMETRY 

There  seems  to  be  no  standard  method  used  in  the 
literature  to  define  ogive  shapes.  An  ogive  shape  is  a 
pointed  body  formed  by  rotating  a circular  arc.  If  the 
water-impact  forces  do  indeed  depend  chiefly  upon  the 
geometry  of  the  body,  the  shape  will  have  to  be  defined 
more  precisely.  The  method  chosen  for  use  in  this  report 
uses  two  angles  and  one  reference  length.  One  angle  is 
the  tip  angle  of  a cone  of  the  same  length  and  base  as 
the  ogive  and  the  other  angle  is  the  arc  length  (swept 
angle)  used  to  generate  the  ogive  section.  The  base 
diameter  is  used  as  the  reference  length.  The  tip  angle, 
base  angle,  slopes,  and  radii  are  all  easily  computed 
for  both  normal  ogives  and  inverse  ogives  (cusps) , as 
shown  in  Figure  A-l.  Note  that  in  Figure  A-l  the  inclusion 
of  a minus  sign  with  the  swept  angle  0 changes  the  ogive 
to  the  cusp.  An  inspection  of  the  computer  program  shows 
a continuing  similarity  in  the  geometric  equations.  The 
only  difference,  if  any,  is  a minus  sign  for  the  cusp 
shape.  The  swept  radius  r is  one  parameter  that  has  the 
same  equation  for  both  shapes. 


r = r/2  sin  (ct/2)  sin  (6/2) 


In  the  comparison  of  the  results  for  different  ogives, 
it  is  often  possible  to  relate  the  forces  to  the  base  area 
of  a tangent  ogive.  A method  for  computing  the  ratio 
between  the  base  area  of  a secant  ogive  and  corresponding 
tangent  ogive  is  shown  in  Figure  A-2.  Also,  the  nomen- 
clature must  be  such  that  the  actual  tip  angle  remains 
constant. 

In  this  report  it  is  assumed  that  the  tangent  ogive  is 
a limit  of  this  geometric  family  of  shapes.  This  implies 
that  the  swept  angle  may  not  be  larger  than  the  cone  angle. 
It  is  also  defined  that  the  limit  case  for  cusps  is  when  the 
swept  angle  is  equal  to  the  cone  angle.  For  cone  angles  of 
90  degrees  or  more,  this  definition  merely  states  that  the 
base  angle  cannot  be  negative. 
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APPENDIX  B 

This  section  consists  of  a listing  and  sample  run  of  the 
computer  program  written  in  BASIC  used  to  calculate  the 
predicted  values  of  the  drag  coefficient  given  in  this  report. 
By  the  selection  of  suitable  constants,  it  is  possible  to 
obtain  predicted  drag  coefficient  versus  depth  functions  for 
shapes  defined  by  the  polynomial 

R = Kl*D  + K2*D2  + K3*D3  + K4*D4  + K5*D5  + K6*D*3 

or  for  ogives  or  cusps.  It  is  also  possible  to  obtain  the 
printout  either  at  the  point  of  the  first  maximum  of  the  total 
drag  coefficient  or  at  uniformly  spaced  distance  intervals. 

The  format  used  in  this  program  was  chosen  in  an  effort  to 
allow  ready  understanding,  easy  modification,  and  maximum 
flexibility. 


The  major  sections  of  the  program  are: 

Lines  100-180  consist  of  an  alphabetical  list  of  constants 
needed  to  operate  the  program. 

Lines  200-550  contain  the  commands  that  establish  the  order 
of  calculations.  Also  the  penetration  ratio  is  calculated  on 
line  250. 

Line  554  contains  the  print  command  used  with  the  condensed 
print  format. 

Lines  700-790  contain  the  calculation  of  the  radius  and  angle 
for  shapes  defined  by  the  polynomial. 

Lines  900-990  contain  the  calculation  of  the  radius  and  angle 
for  ogives  and  cusps.  The  headings  used  with  these  shapes  are 
also  printed. 

Lines  1200-1290  contain  the  calculations  of  the  skin  friction 
drag  coefficient.  A constant  value  of  .003  is  assumed  for  the 
friction  drag  coefficient. 

Lines  1400-1490  contain  the  calculations  of  steady-state 
pressure  drag  coefficient. 

Lines  1600-1690  contain  the  calculations  of  the  transitory 
drag  coefficient. 

Lines  2000-2290  contain  the  calculations  of  transitory 
and  total  drag  coefficients  provided  the  base  of  the  forebody 
shape  is  submerged. 
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Lines  8000-8100  contain  the  print  statements  used  with  the 
long  print  format. 

In  order  to  run  a problem  it  is  necessary  to  enter  the 
following  constants.  Note:  all  distances  must  have  the  same 

units  and  determine  the  unit  used  in  output. 

1.  Line  115  Dl  = distance  between  computational  step. 

In  long  print  format  every  tenth 
computational  step  is  printed 

2.  Lines  135-151  K1-K6  = coefficients  of  polynomial 

Lines  152  K1  = total  cone  angle  in  degrees  of  ogive  or  cusp 
Lines  153  K8  = swept  angle  in  degrees  of  ogive  or  cusp 

3.  Lines  154  P = -1  for  cusp 

P = 1 for  ogive 

P = 2 for  polynomial 

4.  Line  157  P4  = 0 printout  of  first  maximum  only 

P4  = 1 full  printout 

5.  Line  160  R 1 = radius  of  base 

The  various  constants  are  defined  as: 

A.  temporary  name  for  an  angle 

A.  body  angle  at  the  original  water  surface 

A3  previous  A4 

A4  body  angle  at  the  effective  surface 

C C6  when  R4=R1 

Cl  incremental  friction  drag  coefficient 

C2  friction  drag  coefficient 
C3  previous  C4 

C4  steady-state  pressure  drag  coefficient 
C5  previous  C6 

C6  transitory  drag  coefficient 
C7  previous  C8 

C8  total  drag  coefficient  (C2  + C4  + C6) 

Dl  distance  increment 

D2  distance  from  tip  to  original  water  surface 
D3  Dl/H 

D4  distance  from  tip  to  effective  water  surface 

H penetration  ratio 

K1-K6  constants  that  determine  the  polynomially  defined 
shapes 

K7  cone  angle  of  ogive  or  cusp 
K8  swept  angle  of  ogive  or  cusp 
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radius  of  curvature  of  ogive  or  cusp 
base  radius 

body  radius  at  original  water  surface 
previous  P4 

body  radius  at  effective  water  surface 

temporary  area 
base  area 

area  at  ori9inal  water  surface 

surf!™!  Jtl0n-area  at  effective  water  surface 
surface  area  increment 
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LIST  Fig.  B-l  - Listing  and  Sample  Run  of  Computer 

Program  Written  in  BASIC 

75/07/10.  17.19.18. 


PROGRAM 

OBJ  46  4 

1 01 

LET 

A4= . 000 1 

1 02 

LET 

C2  = 0 

1 05 

LET 

C4=0 

1 10 

LET 

C8  = 0 

1 15 

LET 

D1  = • 0 1 

1 20 

LET 

1)2  = 0 

1 25 

LET 

D4=0 

1 30 

LET 

D5=0 

1 34 

LET 

K=  1 

1 35 

LET 

K l =0 

1 40 

LET 

K2=0 

1 45 

LET 

K3=0 

1 49 

LET 

K4=0 

1 50 

LET 

K5= . 38 

1 5 1 

LET 

K6=0 

152 

LET 

K7=60 

153 

LET 

K8=  43 

15  4 

LET 

P=1 

155 

LET 

Pl=-1 

156 

LET 

P3=  1 

15  7 

LET 

P4=2 

1 58 

LET 

P5=0 

1 60 

LET 

Rl  = 1 

1 65 

LET 

R3=0 

1 66 

LET 

R4=0 

1 70 

LET 

Sl=3. 14159*Rlt2 

180 

LET 

Z=0 

2 00 

LET 

P1=P1+1 

2 09 

IF  P4=0  THEN  220 

2 10 

IF  P 1 <9 .5  THEN220 

2 15 

G0  SUB  8000 

220 

LET 

D2=D2+D1 

2 30 

LET 

D=D2 

235 

IF  P< 1 • 5 THEN  245 

2 40 

G0  SUB  700 

2 41 

REM 

SUB  900  IS  FOR  0 

2 43 

G0  T0  250 

2 45 

G0  SUB  900 

2 50 

LET 

H = 1 - . 6905  E"2*A+8 

260 

LET 

A2=  A 

2 70 

LET 

R2=R 

280 

LET 

D3=D1/H 

290 

LET 

D4=D4+D3 

300 

LET 

D=D4 

3 10 

LET 

R3=R4 

320 

LET 

A3=A4 

3 25 

IF  P< 1 • 5 THEN  335 

3 30 

G0  SUB  700 

CONT. 
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3  32  G0  T0  340 

3  35  G0  SUB  900 

3  40  LET  A4=  A 

3  50  LET  R4=R 

3  60  IE  R4>R1  THEN  2000 

3 70  IF  A4<  • 1 THEN  2000 

380  G0  SUB  1200 

385  G0  SUB  1400 

390  IF  A4>A2  THEN  420 

4 00  LET  A=A2 

4  10  G0  T0  45  0 

420  LET  A=A4 

450  G0  SUB  1600 

4  60  IF  P5>.5  THEN  470 

461  LET  C5=C2+C4 

462  IF  C3>C5  THEN  465 

4  63  LET  C3=C 5 

464  G0  T0  470 

4  65  LET  P5=  1 

4 70  LET  C7  = C8 

480  LET  C8=C3+C6 

490  IF  C8<  • 02  THEN  200 

5 00  G0  T0  510 

510  IF  C3<C4  THEN  530 

520  LET  C4=C3 

5  30  LET  C8  = C3+C6 

535  IF  P4=0  THEN  540 

5  36  G0  T0  550 

5  40  IF  C7<C8  THEN  550 

5  44  PRINT  D2/C8/C3/A2/D4 

5  46  G0  T0  9999 

5 50  G0  T0  200 

700  LET  R=Kl*D+K2*Dt2+K3*Dt3+K4*D»4 
7 05  LET  R=R+K5*Dt  . 5+K6*D»  .3 

710  LET  A=Kl+2*K2*D+3*K3*Dt2+4*K4*Dt3 

7 15  LET  A=A+.5*K5*DtC-.5>+.3*K6*DtC-.7> 

720  LET  A=57.3*ATN< A) 

7 60  IF  P3>  1 THEN  790 

7 62  PRINT  ’’SHAPE  GIVEN  BY  POLYNOMIAL” 

763  PRINT  ”P0LY  CONST  ="K1,K2,K3 

764  PRINT  K4,K5#K6 

765  PRINT 

766  PRINT  ”B0DY  RADIUS=”R1 

767  PRINT 

768  PRINT  ”UIST”,”T0TAL  CD”/  ’’CDS”# ’’ANGLE  ”/ ”EFF  DIST” 

785  LET  P3=2 

790  RETURN 

9 00  IF  1 THEN  930 

9 05  LET  U=  • 5*K7/57 • 3 

9 10  LET  V=.5*K8/57.  3 

920  LET  Q=R1/<C0S<U-V)-C0S<U+V>> 

9 30  LET  X = 0*SIN(  U+P*V  > -P*D 


CONT. 
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940  LET  Y=SQR(Gt2-Xt2) 

9 50  LET  A=ATN( X/Y> 

9 60  LET  R=P*Q*(C0SCA)-C0S(  U+P*\/  ) ) 

9 70  LET  A=57.3*A 
9 71  F P3>  1 THEN  99  0 
9 72  IF  P>-. 5 THEN  975 
9 73  PRI NT"THE  SHAPE  IS  A CUSP'* 

9 74  G0  T0  976 

9 75  PRINT  "THE  SHAPE  IS  AN  OGIVE" 

9 76  PRINT 

9 77  PRINT  "BODY  RADIUS="R1 
9 78  PRINT  "C0NE  ANGLE  ="K7 

979  PRINT  "SWEPT  ANGLE  ="K8 

980  PRINT  "RAD  0F  CURV  ="Q 

981  PRINT 

982  PRINT"DIST'S  "TOTAL  CD'S  "CDS'S  "ANGLE'S  "EFF  DIST" 

988  LET  P3=2 

990  RETURN 

1 200  LET  S6=3.  141  59*<  R3+R4)  *D3 
1210  LET  C 1 = » 003+S6/S 1 
1 220  LET  C2=C2+C 1 
1290  RETURN 
1 400  LET  C3=C4 

1 410  L-.  T C4=3E-3- 1 • 36E-3*A+ 1 • 076E-3*At 2-3»04E-5*At3+3. 41 5F- 7*a»4 
1420  LET  C4=C4-1.362E-9*A»5 

1430  LET  S=3. 14159*R2»2 
1440  LET  C4=C4*S/S1 
1450  LET  C4=C4/H»2 
1 490  LET  C9=C4 

1600  LET  C6=  1 • 76E-2+1  • 093E-2*A-8  • 663E- 4*A»  2+3*9  12E -5*A»3 
1610  LET  C6=C6-2.556E-7*A»4 
1620  LET  C6=C6*TAN((A-3)/57.3) 

1630  LET  C6=C6*S/S1 
1690  RETURN 

2 000  PRINT 

2005  LET  D2=D2- D1 
2010  G0  SUB  8000 
2015  IF  P4<  • 5 THEN  9000 
2 020  LET  A=A2 
2 030  IF  P>  • 5 THEN  2050 
2 040  LET  A=A4 
2 05  0 LET  R=R1 
2 060  LET  C=C6 
2 065  LET  D=D2 
2 067  LET  T9  = A/57.3 
2 068  LET  T8  = TAN(  T9) 

2 070  LET  M2=D2 
2 080  LET  T2=10*D1 
2 200  LET  B4=  •661  *T2/M2 

2210  LET  A9  = »3/<  1 + C T2/R I ) ) 1 2+EXPC  - 9*B4) 

2220  LET  S6=C*(  1/1.3)*A9 

2230  LET  C8=C3+S6 

2240  UCT  D2=T2+D 
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2 250  PRINT  D2> C8 

2260  LET  T2=T2+10*D1 

2262  IF  T2<  »8*D  THEN  2270 

2 265  LET  T2=T2+50*D1 

2270  IF  S6< « 02*C3  THEN  9000 

2 290  G0  T0  2200 

2990  G0  T0  9000 

8 00u  PRINT  D2>C8#C3>  A2>  D4 

8 050  LET  P1=0 

8060  LET  Z = Z+1 

8 065  IF  Z<4«  5 THEN  8 100 

8070  PRINT 

8 075  LET  Z=0 

8 100  RETURN 

9 000  LET  Z = Z 
9999  fMD 
READY. 


LOG  0FF. 
CP 


17.26.52. 
1.518  SEC. 


k' 
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RUN 


75/07/10.  17.16.26. 
PROGRAM  OBJ  46  4 


THE  SHAPE  IS  AN  0GIV E 

BODY  RADI  US=  1 
CONE  ANGLE  = 60 
SWEPT  ANGLE  = 43 
RAD  OF  CURV  = 2.72888 


DIST 

. 1 
.2 
.3 
. 4 
.5 

* 6 
.7 
.8 
.9 
1 . 

1 . 1 
1.2 

1.3 

1 . 4 

1 . 47 

1 . 57 
1.67 
1 . 77 
1 .87 
1.97 


CP  1.099 


TOTAL  CD 
3. 33095E-2 
9. 71 428E-2 
• 1 6 1 6 28 
•214612 
.252057 

• 2736  42 

• 280685 
•275171 
•261225 
•247124 

•23441 7 
•223612 

• 2 1 48  1 
•207832 

• 2038  52 

• 198672 

• 195  15  1 

• 192743 

• 1 9 1 08  4 
•189931 


SECS. 


CDS 

1 . 15 109E-2 
•038293 
7.1  7684E-2 

• 106033 

• 136926 


ANGLE 
48.2419 
AS.  1 79  4 
42.2736 
39.4961 
36.8255 


• 14  1 58  1 
. 1 78  18  7 

• 185841 

• 186357 

• 18635  7 


34.245 

31.7414 

29.3037 

26.9229 

24.5914 


. 18635  7 

• 186357 
. 186357 

• 186357 


22.3025 
20.0506 
1 7.8305 
15.6378 


•186357  13.9005 


E FF  DIST 
•126227 
.250941 

• 37427 
•496305 
. 6 1 7 1 08 

.736721 

• 8551  73 

• 9 72478 
1 .08864 
1 .20366 

1.31 754 
1.43025 
1 .541  78 
1 .6521 

1 .739  49 


RUN  COMPLETE 


'■  -’iff » 


APPENDIX  C 

This  section  consists  of  graphs  which  compare  the 
experimental  and  predicted  values  of  total  drag  coefficient 
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PREDICTED 

EXPERIMENTAL  REF  fe) 


. C-2  TOTAL  DRAG  COEFFICIENT  VS  DEPTH  FOR  60/60  OGIVES 
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PREDICTED 


FIG.  C-4  TOTAL  DRAG  COEFFICIENT  VS  DEPTH  FOR  41 .44/41 .44  OGIVES 
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